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Abstract 

We suggest that galileon theories should have an additional self-coupling of the fields to the trace of their own energy- 
momentum tensor. We explore the classical features of one such model, in flat 4D spacetime, with emphasis on solutions that 
are scalar analogues of gravitational geons. We discuss the stability of these scalar geons, and some of their possible signatures, 
including shock fronts. 
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Galileon theories are a class of models for new scalar 
fields whose Lagrangians involve multilincars of first and 
second derivatives, but whose nonlinear field equations 
are nonetheless still only second order. They may be 
important for the description of large-scale features in 
astrophysics as well as for elementary particle theory [l|, 
7] . Hierarchies of such Lagrangians giving rise to such 
field equations were first discussed mathematically in 0- 
10, 14]. The simplest example involves a single scalar 
field. 

This galileon field is usually coupled to all other mat- 
ter through the trace of the energy-momentum tensor, 
q (matter)^ anc j j g tf^g gravitation-like by virtue of the 
similarity between this universal coupling and that of 
the metric g M „ to 0^™ attGI ' ) in general relativity. Indeed, 
some galileon models have been obtained from limits of 
higher dimensional gravitation theories Q . 

But surely, in a self-consistent theory, for the coupling 
to be truly universal, the galileon should also be cou- 
pled to its own energy-momentum trace, even in the flat 
spacetime limit. Some consequences of this additional 
self-coupling are considered in this paper. 

The action for the lowest non-trivial member of the 
galileon hierarchy can be written in various ways upon 
integrating by parts. Perhaps the most compact and 
memorable of these is 



A 2 = \ 



x . 



(1) 



where <j> is the scalar galileon field, 4> a = dcf> (x) / dx a , 
etc., and where repeated indices are summed using the 
Lorentz metric 5^ v = diag (1, — 1, — 1, ■ ■ ■ ). 

It is straightforward to include in A2 & covariant cou- 
pling to a background spacetime metric and hence to 
deduce a symmetric energy-momentum tensor. In the 
flat-space limit, the result is 

= <t>ii<t>^4>aa - <l>a<t>av<t>n ~ 4>a4>a^4>v + Spv^u^p^aP ■ 

This is seen to be conserved, 



8 (2) 



(3) 



upon using the field equation that follows from locally 
extremizing A%, = SA2/5(f> — — £ 2 [</>], where 



(2) 

An interesting wrinkle now appears: 0^ is not trace- 
less. Consequently, the usual form of the scale current, 
x a Oal, i s n °t conserved [l5[. On the other hand, the 
action |T]) is homogeneous in <f> and its derivatives, and is 
clearly invariant under the scale transformations x — > sx 
and 4>{x) — > s( 4 ~"'/ 3 (sx). Hence the corresponding 
Noether current must be conserved. This current is eas- 
ily found, especially for n = 4, so let us restrict our at- 
tention to four spacetime dimensions in the following. 

In that case the trace is obviously a total divergence: 



9 (2) 
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9 a {4>a4>p4>/3) 



(5) 



That is to say, for n — 4 the virial is the trilinear V a = 
4>a4>p4>P- So a conserved scale current is given by the 
combination, 



(6) 



Interestingly, this virial is not a divergence modulo a con- 
served current, so this model is not conformally invariant 
despite being scale invariant. Be that as it may, it is not 
our principal concern here. 

Our interest here is that the nonzero trace suggests an 
additional interaction where <f> couples directly to its own 
O*- 2 ). This is similar to coupling a conventional massive 
scalar to the trace of its own energy-momentum tensor 
[Til ] . In that previously considered example, however, 
the consistent coupling of the field to its trace required 
an iteration to all orders in the coupling. Upon sum- 
ming the iteration and making a field redefinition, the 
Nambu-Goldstonc model emerged. But, for the sim- 
plest galileon model in four spacetime dimensions, ([T]), a 
consistent coupling of field and trace is much easier to 
implement. No iteration is required. The first-order 
coupling alone is consistent, after integratin g by parts 
and ignoring boundary contributions, so that [161 ] 



(f> d a (cj) a (l)p(f)p) d 4 x 



x . 



(7) 



£2 [<t>] = fiaa&PP - fyap^aP 



(4) 



(Similar quadrilinear terms have appeared previously in 
d, |1] , only multiplied there by scalar curvature R so 
that they would drop out in the flat spacetime limit that 
we consider.) Consistency follows because ([7]) gives an 
additional contribution to the energy-momentum tensor 
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which is traceless, in 4D spacetime: 



©ft 3 !/ = 4>n^v4>a^ot 



\&liv<Pa<Pa. 



= o . (8) 



Of course, coupling <p to its own trace may impact the 
Vainstein mechanism [2(j by changing the effective cou- 
pling of ©( mattor ) to both backgrounds and fluctuations 
in (p. We leave this as an exercise for the reader. 

Based on these elementary observations, we consider a 
model with action 

A = J {\(t>a<t><x - \Ma<$>a4>PP ~ \^<Pa4'a<p0 ( Pp) ^X , 

where for the Lagrangian L we take a mixture of three 
terms: the standard bilinear, the trilinear galileon, 
and its corresponding quadrilinear trace-coupling. The 
quadrilinear is reminiscent of the Skyrme term in non- 
linear a models [l9| although here the topology would 
appear to be always trivial. 

The second and third terms in A are logically con- 
nected, as we have indicated. But why include in A 
the standard bilinear term? The reasons for including 
this term are to soften the behavior of solutions at large 
distances, as will be evident below, and also to satisfy 
Derrick's criterion for classical stability under the rescal- 
ing of x. Without the bilinear term in L the energy 
within a spatial volume would be neutrally stable under 
a uniform rescaling of x, and therefore able to disperse 

Similarly, for positive k, the last term in A ensures the 
energy density of static solutions is always bounded below 
under a rescaling of the field (p, a feature that would not 
be true if K = but A =/= 0. So, we only consider k > in 
the following. But before discussing the complete 
for the model, we note that we did not include in A a 
term coupling (p to the trace of the energy-momentum 
due to the standard bilinear term, namely, J (pQ^d 4 x, 
where 



Qui = §\i§v - h&nv<Pu<Pa , (1) = -<t>a 



(10) 



We have omitted such an additional term in A solely as 
a matter of taste, thereby ensuring that L is invariant 
under constant shifts of the field. Among other things, 
this greatly simplifies the task of finding solutions to the 
equations of motion. 

The field equation of motion for the model is = 
5 A/ 5(f) = —£ [</>], where 

£ [</>] = (j) aa - A (4> aa (j)i3p - <p a $(pa(i) - K (<j) a 4>P<l>p) a ■ 

(11) 

As expected, this field equation is second-order, albeit 
nonlinear. Also note, under a rescaling of both x and 
(p, nonzero parameters A and k can be scaled out of the 
equation. Define 



(12) 



Then the field equation for ip (z) becomes 

VW - (^aa-0/3/3 ~ Ipafitpap) ~ (V'a##) Q = , (13) 

where ip a = dip(z)/dz a , etc. In effect then, if both 
A and n do not vanish, it is only necessary to solve the 
model's field equation for A = K = 1. 

Though £ is nonlinear, it is nevertheless still true that 
some plane waves are exact solutions. For "light-ray" 
plane waves, £ [Aexp (ik a x a )] = for constant A and 
k a , if k a k a = with A arbitrary. In this case, each 
of the terms in £ vanish separately. In fact, light-ray 
plane waves are only one among many possible solutions 
for which both 4> aa = and tf>p<t>p — 0. On the other 
hand, for massive plane waves, £ [Aexp (ik a x a )] = if 
l/k a k a = —3kA 2 < 0. The latter "tachyonic" solutions 
would seem to be less interesting for real physics. 

For static, spherically symmetric solutions, <fi — 4>{r), 
the field equation of motion becomes 



u 1 d 

r 2 dr 



4>' + \-{cj> l f + k{4>')- 



(14) 



where <j>' = dip /dr. This is immediately integrated once 
to obtain a cubic equation, 



rV + 2Ar (<f/f + nr 2 (ip'f = C , 



(15) 



where C is the constant of integration. Now, without 
loss of generality (cf. (fl2|) and (1X31 ) we may choose A > 0. 
Then, if C = 0, either <p' vanishes, or else there are two 
solutions that are real only within a finite sphere of radius 
r = yj\ 2 /k. These two "interior" solutions are given 
exactly by 



rn V 



(16) 



Note that these solutions always have ip' < within the 
finite sphere. 

Otherwise, if C ^ 0, then examination of the cubic 
equation for small and large \<f/\ determines the asymp- 
totic behavior of <p' for large and small r. In particular, 
there is only one type of asymptotic behavior for large r: 



C 

(pi ~ -^r for either sign of C 

r—>oo t 

However, there are two types of behavior for large 
corresponding to small r. Either 

-2A 



4>'k 



i + o[1 



(17) 



(18) 



provided (p 1 < 0, but with either sign of C; or else 



(19) 



provided C > 0, but with either sign of <p' . The corre- 
sponding real solutions behave as 

— 2A 

for either sign of C, or (20) 

(21) 



i"->0 KT 



C 



±\l — — provided C > 
r^o V 2Ar 1 
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Comparison of the small r behavior to the large r asymp- 
totics shows that in half these cases the solutions would 
require zeroes to be real and continuous for all r. But 
such zeroes do not occur. Instead, half of the cases pro- 
vide real solutions only over a finite interval of r, some- 
what similar to the C — solutions in (fl"6|) . but not so 
easily expressed, analytically. 

The solutions which are real for all r > boil down 
to two cases, with small and large r behavior given by 
either 

0' ~ \ WT- and 4>' ~ 4 forC*>0, (22) 
i — >o V 2Ar ' — s-oo r 

or else 

— 2A C 
(/)' and 0' ~ — forC<0. (23) 

From further inspection of the cubic equation to deter- 
mine the behavior of <\>' for intermediate values of r, when 
C > it turns out that <fi' is a single- valued, positive 
function for all r > 0, joining smoothly with the asymp- 
totic behaviors given in (|22|) . However, it also turns 
out there is an additional complication when C < 0. 
In this case there is a critical value (k 3//2 /A 2 ) C cr iti ca i = 
-4^/27 fa -0.2566 such that, if C < C critica i then </>' 
is a single- valued, negative function for all r > 0, while 
if Ccriticai < C < then 0' is triple-valued for an open 
interval in r > 0. It is not completely clear to us what 
physics underlies this multivalued-ness for some negative 
C. But in any case, when C < it is also true that 
<fJ joins smoothly with the asymptotic behaviors given 
in (|23|) . All this is illustrated in Figures 1 and 2, for 

A = K = 1. 

A test particle coupled by <^Q( matter ) to any of these 
galileon field configurations would see an effective po- 
tential which is not 1/r, for intermediate and small r. 
Therefore its orbit would show deviations from the usual 
Kepler laws, including precession at variance with that 
predicted by conventional general relativity. It would be 
interesting to search for such effects, say, by considering 
stars orbiting around the galactic center. 

For the solutions described by (|2"21 and (|2"3"j) . the to- 
tal energy outside any large radius is obviously finite for 
both C > and C < 0. And if C > 0, the total energy 
within a small sphere surrounding the origin is also man- 
ifestly finite. But if C < the energy within that same 
small sphere could be infinite unless there is a cancella- 
tion between the galileon term and the trace interaction 
term. Remarkably, this cancellation does occur [l7| . So 
both C > and C < types of static solutions for the 
model have finite total energy. 

Complete information about the distribution of energy 
is provided by the model's energy-momentum tensor, 

= - Ae; t 2 J - «e$ . (24) 

As expected, this is conserved, given the field equation 
£[</>] =0, since 

B„Q lu/ = <f> v ey>] . (25) 
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FIG. 1: tp'(r) for C = +1/4^, with AT = 0,1, 2, 3 for top to 
bottom curves, respectively. 




FIG. 2: ip' (r) for C = -1/2 N , with N = 6,5, 4, 3, 2, 1, from 
left to right, respectively. The thin black curve is a union of 
the two C — solutions in (|16[) . 



The energy density for static solutions differs from the 
canonical energy density for such solutions (namely, —L) 
by a total spatial divergence that arises from the galileon 
term: 

©00 = - £| static ~ ^ ■ ((V0) 2 ^) . (26) 

This divergence will not contribute to the total energy for 
fields such that linv^oo (<f>/\a.r) exists. Assuming that 
is the case, Derrick's scaling argument for static, finite 
energy solutions of the equations of motion [4| shows the 

energy is just twice that due to the bilinear Qqq. Thus, 

E = J 6 00 d 3 r = J (V</>) 2 d 3 r . (27) 

For the spherically symmetric static solutions of (|15l) . 
this becomes an expression of the energy as a function of 
the parameters and the constant of integration C: 

/>CO 

E [A, k, C\ = 4tt / (<f)' f r 2 dr . (28) 
Jo 
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E(±C) 
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FIG. 3: E (±C) versus C > as lower/upper curves (the 
horizontal line is E (Critical) ~ 3.7396). 



Again without loss of generality, consider A = k = 1. 
Then for either C > or for C < Critical < [3, 
change integration variables from r to s = </>' to find: 



£(C^0)=/(|C|)t(|C| + ±7t) 



Z(<7> 0) 



P(s,C) (is 
o (s 2 + l) 4 R (s, C) 



(29) 
(30) 



where R (s, C) — ^ s 4 + s (s 2 + 1) C and where the nu- 
merator of the integrand is an eighth-order polynomial 
in s, namely, P (s, C) = 8s 8 + YlCs 1 + (3C 2 - 8) s 6 + 
8Cs 5 + 7C 2 s 4 - 4Cs 3 + 5C 2 s 2 + C 2 . Thus, I (C) is an 
elliptic integral. But rather than express the final result 
in terms of standard functions, it suffices here just to plot 
E (C), in Figure 3. Note that E increases monotonically 
with \C\. 

For other values of A and n with the constant of inte- 
gration C specified as in (Ti~5|) , the energy of the solution 
is given in terms of the function defined by (|29|30p : 

E[X,k,C] = (a 3 /k 5/2 ) e(k 3/2 C/A 2 ) . (31) 

The energy curves indicate double degeneracy in E, for 
different values of |C|, when E[X,k,C] > 7tA 3 /k 5 / 2 . 
Also, for a given |C| the negative C solutions are higher 
in energy, with E [A, re, - |C|] - E [A, k, \C\] = ttA 3 /k 5 / 2 + 
2 |C| A/k. Or at least this is true for all \C\ > |C cr iti ca i| 



> 



,E 



-a 



critical 



in which case E [A, K, C] 

3.7396 A 3 /k 5 / 2 

Finite energy classical solutions of gravity-like theo- 
ries brin g to mind the "geons" proposed long ago by 
Wheeler [2l[. These were envisioned in their purest 
form as distributions of only gravitational energy held to- 
gether solely by gravitational interaction. Combinations 
of electromagnetic energy and gravity were also consid- 
ered, as were systems containing neutrinos. Wheeler 
argued that such configurations would be relatively sta- 
ble, if they existed, but would eventually dissipate due 



to a variety of both classical and quantum effects, in- 
cluding light-light scattering, as well as production and 
absorption of quanta. While plausible distributions were 
sketched, and decay rates were estimated, exact classical 
solutions were not found. 

The same mechanisms would seem to apply to any hy- 
pothetical classical galileon distributions such as those 
discussed here, the main difference being that analytic 
spherically symmetric solutions might still be obtainable 
_ even if conventional gravitational effects were included. 

Perhaps these gravitational effects would not alter the 
, qualitative features of the static pure <f> configurations 
5 given above. Should they really exist, presumably these 
galileon geons could also be dissipated by various classical 
and quantum effects. All this is far beyond our current 
abilities and the scope of this paper, of course, but the 
general ideas suggest some interesting possibilities. 

Whatever the cause, if the configuration's energy loss 
were gradual, as a first step it might suffice to model 
the time-dependent system quasi-statically, as a contin- 
uous flow from one static solution to another. That is 
to say, perhaps a good approximation would be to take 
C (t), with |C| and E (C) decreasing monotonically with 
time. For the positive C case, this would be more or less 
uneventful as the whole configuration would just slowly 
disappear without any abrupt changes. But for the neg- 
ative C CcLSG, clS t increased Critical would be reached, 
beyond which the solution would begin to fold over, ex- 
hibiting the multivalued features shown in Figure 2. But 
this is just the usual picture for the formation of a shock 
front. These particular galileon shocks would implode, 
converging towards the origin, as shown here. We be- 
lieve this is a plausible scenario and a reasonable phys- 
ical interpretation of the model's multivalued solutions. 
Moreover, it would seem to provide a signature for their 
existence. 

As is clear from Figure 2, the shock front would form 
when dip' /dr = oo. For the C < static solutions of (JT5J) 
it is not difficult to determine the locus of such singular 
points. It is given by the intersection of the solutions, 



for various C, and the curve (1 + 3kc 



= 4A0'. As 



usual for singular points in the development of a shock, 
almost certainly there is some physics missing from the 
equations. Since 4>" is large, the obvious modification 
would be to include higher derivative terms in the action, 
which is tantamount to attempting an ultraviolet com- 
pletion of the model. This is an open question. Perhaps 
higher terms in the galileon hierarchy would be natural 
candidates to be included. 

To get a handle on such terms, and for purposes of 
comparison to the model in © , consider briefly another 
model somewhat similar in form, but whose Lagrangian 
consists only of terms taken from the galileon hierarchy, 
without any coupling to O. After rescaling the field and 
coordinates to achieve a standard form, this alternate 
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model may be defined by 

^self-dual W\= J (\lpa4>a- \lp a4> a 1p f3 P 

+j^ip a ip a (ipppip^ - ipp^ipp^)) d 4 x . (32) 

The difference with © lies in the last term, which is 
quadrilinear in the field, as before, but now has two fields 
with second derivatives. 

As the name suggests, this model is self-dual, in the 
following sense: The action retains its form under a 
Legendre transformation @ (also see [13) to a new field 
VP and new coordinates X, as defined by: 

ip (x) + * (X) = x a X a . (33) 

Thus A se if. dua i [ip] = A scU -duai [*], provided integrations 
by parts give no surface contributions. This identity sug- 
gests that there are interesting properties for the quan- 
tized model, such as its ultraviolet behavior, but that is 
outside the scope of the present discussion. 

Here it suffices to compare the classical physics follow- 
ing from (f3"2"| with that following from Upon inte- 
grating once the classical equations of motion for static, 
spherically symmetric solutions of the field equations for 
(|3"2")l . the result is again a cubic equation, 

rV + r(^) 2 + jW 3 = C, (34) 

but the (ip') 3 term is no longer weighted by r 2 as it was 
in (|15[) . Thus the small and large r behaviors are now 
given by 

ip' ~ (3C) 1/3 and ip' ~ (35) 

r— >0 r— >oo T 



for either sign of the constant of integration, C. These 
static solutions have finite total energy for either sign of 
C, as before, only now ip' is always bounded. Moreover, 
upon inspection of the behavior of ip' for intermediate r, 
and various C, unlike the previous model the solutions 
are now always single- valued for either C > or C < 0. 
Thus there are no multivalued solutions like those shown 
in Figure 2. However, each of the C < static solutions 
does have a single point for which dip' /dr = oo, namely, 

1 /3 

r = (3 |C|) ' . So there is still a reason to expect the 
existence of shock fronts for quasi-static time-dependent 
fields in this alternate model. Finally, again for C < 0, 
to have <p' real for all r > 0, it is necessary to join together 

1 /3 

"interior" and "exterior" solutions at r = (3|C|/2) ' . 

It remains to investigate the stability of these spheri- 
cally symmetric solutions under perturbations, especially 
to check for the existence of superluminal modes, along 
the lines of [13j. Evidently, superluminal modes are a 
possible feature for models of this type. 

In conclusion, it would be interesting to search for evi- 
dence of geons containing galileons at all distance scales, 
including galactic and sub-galactic, as well as cosmolog- 
ical. Perhaps a combination of trace couplings and var- 
ious galileon terms, such as those in ((9]) and (|32|) . will 
ultimately lead to a realistic physical model. 
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